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TECHNICAL  NOTE  2129 

METHOD  OF  CALCULATING  THE  LATERAL  MOTIONS  OF  AIRCRAFT 
BASED  ON  THE  LAPLACE  TRANSFORM 
By  Harry  E.  Murray  and  Frederick  C.  Grant 

SUMMARY 


The  lateral  motions  of  aircraft  are  obtained  by  means  of  the 
Laplace  transform  which  gives  solutions  expressed  in  terms  of  elementary 
functions  for  the  free  and  forced  motions.  These  equations  permit  the 
calculation  of  the  free  motion  of  an  aircraft  following  any  initial 
condition  or  the  forced  motion  following  the  application  of  constant 
external  forces  and  moments.  These  forced  motions  can  be  used  to  obtain 
by  means  of  Duhamel's  integral  the  response  to  any  arbitrary  forcing 
function.  All  the  classical  stability  concepts  can  be  deduced  from 
these  same  solution  equations  largely  by  inspection.  These  equations 
for  the  lateral  motion  are  applied  to  the  calculation  of  the  lateral 
stability  of  a  specific  airplane  and  to  the  calculation  of  certain  of 
its  free  and  forced  motions. 


INTRODUCTION 


The  lateral  motions  of  aircraft  are  represented  by  three  simultaneous 
differential  equations  which  are  generally  assumed  to  be  linear.  The 
fundamental  problem  of  lateral  dynamics  involves  the  solution  of  these 
differential  equations  in  terms  of  the  aerodynamic  and  mass  parameters 
of  the  airplane.  The  solutions  can  then  be  used  to  obtain  numerically 
the  motion  of  the  airplane  as  a  function  of  time. 

The  recent  application  of  the  Laplace  transform  to  the  solution  of 
systems  of  linear  differential  equations  permits  a  more  general  analysis 
of  the  problem  of  airplane  motion  than  that  of  reference  1,  which  is 
based  upon  Heaviside's  operational  calculus.  Heaviside's  operational 
calculus  permits  a  calculation  of  the  forced  motion,  which  is  the  motion 
following  the  application  of  external  forces  and  moments.  The  Laplace 
transform  permits  these  same  calculations  and  also  permits  the  direct 
calculation  of  the  free  motion,  which  is  the  motion  following  finite 
initial  values  of  the  variables  and  their  first  derivatives  in  the 
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absence  of  externally  applied  forces  and  moments.  This  calculation 
cannot  be  made  by  use  of  Heaviside's  operational  calculus.  The  Laplace 
transform  solutions,  which  include  both  the  free  and  forced  motions,  may 
be  written  in  a  closed  form  from  which  all  the  classical  stability  con¬ 
cepts  can  be  deduced  largely  by  inspection.  The  form  of  the  equations 
of  motions  of  the  airplane  is  independent  of  such  aerodynamic  parameters 
as  Reynolds  number  and  Mach  number,  and  these  parameters  enter  the  equa¬ 
tions  only  as  they  effect  the  values  of  the  aerodynamic  constants  or 
stability  derivatives  appearing  in  the  equations.  The  values  of  the 
stability  derivatives  must  be  obtained  by  actual  measurements  during 
physical  tests  or  from  aerodynamic  theory  before  motion  calculations  can 
be  attempted . 

Investigations  of  some  of  the  possibilities  of  applying  the  Laplace 
transform  to  the  study  of  aircraft  motion  have  been  reported  in  refer¬ 
ences  2  and  3,  and  in  two  British  reports,  one  by  K.  Mitchell,  the  other 
by  J.  Watham  and  E.  Priestley.  The  British  papers  do  not  give  final 
equations  in  a  form  suitable  fbr  calculation  purposes.  The  analysis  of 
reference  2  closely  parallels  that  of  the  present  paper  until  the  point 
of  taking  the  inverse  Laplace  transform  is  reached.  At  this  point, 
reference  2  indicates  that  the  inverse  Laplace  transform  can  be  taken 
either  by  means  of  the  relatively  simple  partial-fraction  expansion 
(used  in  the  present  paper)  or  the  more  complicated  inversion  theorem 
of  the  Laplace  transform.  Neither  approach  in  reference  2  is  carried 
to  the  point  of  final  equations  containing  only  elementary  functions  and 
in  a  form  particularly  suited  for  computation.  A  solution  similar  to 
that  of  the  present  paper  is  indicated  in  reference  3-  Only  the  form 
of  the  analysis  is  shown  in  reference  3,  however,  and  all  the  details 
necessary  for  practical  applications  have  not  been  carried  out. 


The  present  paper  presents  an  analysis  based  on  the  representation 
of  the  lateral  motion  of  an  aircraft  by  differential  equations.  The 
results  of  the  analysis  are  solutions  in  closed  form  expressing  the  free 
and  forced  motions  in  terms  of  elementary  functions.  These  equations 
permit  the  calculation  of  the  free  motion  of  an  aircraft  following  any 
initial  condition  or  the  forced  motion  following  the  application  of 
constant  external  forces  and  moments.  These  forced  motions  can  be  used 
to  obtain,  by  means  of  Duhamel's  integral,  the  response  to  any  arbitrary 
forcing  function  as  shown  in  references  4  and  5-  The  solutions  are 
readily  adaptable  to  calculation  by  digital-type  calculating  machines 
and  the  calculation  is  an  arithmetical  process  requiring  no  knowledge  of 
the  theory  of  the  Laplace  transform.  The  solution  equations  of  motion 
have  been  applied  on  an  automatic  calculating  machine  to  the  calculation 
of  the  lateral  stability  of  a  specific  airplane  and  to  the  calculation 
of  certain  free  and  forced  motions  as  illustrative  examples. 
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COEFFICIENTS  AND  SYMBOLS 

trim  lift  coefficient  (W  cos  7/qS) 

rolling -moment  coefficient  (L/qSb) 

yawing-moment  coefficient  (N/qSb) 

lateral-force  coefficient  (Y/qS) 

airplane  weight,  pounds 

rolling  moment 

pitching  moment 
% 

yawing  moment 
lateral  force 
aileron  hinge  moment 
elevator  hinge  moment 
rudder  hinge  moment 
dynamic  pressure  (pV2/2) 
wing  area,  square  feet 
wing  span,  feet 

inclination  of  flight  path  to,  horizontal  (positive 
in  climh),  degrees 

angle  of  attack,  degrees 

angle  of  pitch,  degrees 

mass  density  of  air,  slugs  per  cubic  foot 

free-stream  velocity,  feet  per  second 

airplane  mass,  slugs  (w/g) 

acceleration  due  to  gravity,  feet  per  second  per 
second 


g 
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sb 

t 


Mb 

♦ 

P 


r 

P 


cz 


c 


C 


nc 


°YC 

P,P' 

Tl/2^Tl/2' 


Rl/2>Nl/2 

6r 


nondimensional  time  (tv/b) 
time,  seconds 


inclination  of  principal  longitudinal  axis  of  inertia 
(positive  for  axis  above  flight  path  at  nose), 
degrees 

airplane  relative-density  factor  (m/pSb) 


angle  of  bank,  radians 


angle  of  yaw  or  azimuth,  radians  (f  r  dty 

rolling  velocity  about  stability  X-axis,  radians  per 
second 

yawing  velocity  about  stability  Z-axis,  radians  per 
second 

angle  of  sideslip,  radians 

rolling-moment  coefficient  of  forcing-function  couple 
in  roll 

yawing -moment  coefficient  of  forcing-function  couple 
in  yaw 

lateral-force  coefficient  of  lateral  forcing  function 

periods  of  oscillatory  modes,  seconds 

times  to  damp  to  half -amplitude  of  oscillatory  modes, 
seconds 

cycles  to  damp  to  half -amplitude  of  oscillatory  modes 


aileron  deflection,  degrees 
rudder  deflection,  degrees 
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Be 

*X 


*Z 


KXZ 


% 


CT 


A 

^"1  *  ^2  >  ^3  >  ^4 


po,qo 

o'  ,  12s 
40  do 


*n 


elevator  deflection,  degrees 

nondimens ional  radius  of  gyration  about  stability 
X-axis 


^KX02cos2n  +  Kz02Bin2ri 


nondimensional  radius  of  gyration  about  stability 
Z-axis 


n  ^cos^n  +  K 
Z0 


*0 


sin 


nondimensional  product  of  inertia  between  stability 
X-  and  Z-axes 


^KZ02  ‘  KX02)  sin  T)  cos  r 


nondimensional  radius  of  gyration  about  principal 
X-axis  (kxo/b) 

nondimensional  radius  of  gyration  about  principal 
Z-axis  (kzQ/b) 


radius  of  gyration  about  principal  X-axis,  feet 
radius  of  gyration  about  principal  Z-axis,  feet 
Laplace  transform  of 
stability  quartic 
roots  of  A  =  0 


polynomials  in  0 


roots  of  q<J  =  0 


real  part  of  X^  and 
complex  conjugates 

x.4 

when 

X3 

and 

are 

imaginary  part  of  X3  and 
are  complex  conjugates 

X^ 

when 

X3 

and 

real  part  of  and 

complex  conjugates 

^2 

when 

X± 

and 

^2 

are 
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imaginary  part  of  an(^-  ^2  w^en  an<^  ^2 

are  complex  conjugates 


A,B,C,D,E 


coefficients  of  stability  quartic 


R  Routh!s  discriminant 

AlfA2,k^,k^,A^,k^  amplitude  coefficients  for  0 


Bi,B2,B3,B4,B5,Bg  amplitude  coefficients  for  \|r 
Cl,C2,C3,C4,C5  amplitude  coefficients  for  3 


RA 

JA 


IB 


RC 


*c 


V 


V 


rb' 


V 


real  part  of  A3  and  A4  when  A3  and  A4  are 

complex  conjugates 

imaginary  part  of  A3  and  A  4  when  A3  and  A4 

are  complex  conjugates 

real  part  of  B3  and  B4  when  A3  and  A4  are 

complex  conjugates 

imaginary  part  of  B3  and  B4  when  A3  and  A4 

are  complex  conjugates 

real  part  of  C3  and  C4  when  A3  and  A4  are 

complex  conjugates 

imaginary  part  of  C3  and  C4  when  A3  and  A4 

are  complex  conjugates 

real  part  of  Aj_  and  A2  when  A-^  and  A^  are 

complex  conjugates 

imaginary  part  of  A]_  and  A2  when  A-j_  .  and  A2 

are  complex  conjugates 

real  part  of  B]_  and  B2  when  A^  and  A2  are 

complex  conjugates 

imaginary  part  of  B^  and  B2  when  A^  and  A2 

are  complex  conjugates 
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V 

V 


kb 


Kc 


ka' 


♦ 


kb' 


V 


cdb 


real  part  of  Cp  and  C2  when  Xp  and  X2  are 
complex  conjugates 

imaginary  part  of  Cp  and  C2  when  Xp  and  X2 
are  complex  conjugates 

amplitude  coefficient  for  0  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  X3  and  X4 

(2  \)rA2  +  Ia2) 

amplitude  coefficient  for  i|r  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  X3  and  X4 

(2  \/Rb2  +  Ib2) 

amplitude  coefficient  for  3  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  X3  and  X4 

(2  vfec2  +  ic2) 

amplitude  coefficient  for  0  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  Xp  and  Xg 

(aA'®  + V2) 

amplitude  coefficient  for  t  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  Xp  and  X2 

(2^-2  +  IB'2) 

amplitude  coefficient  for  3  oscillation  corre¬ 
sponding  to  complex  conjugate  roots  Xp  and  X2 

(2/^-2  +  IC'2) 

phase  angle  for  0  oscillation  corresponding  to 
conjugate  complex  roots  X3  and  X4,  radians 


phase  angle  for  oscillation  corresponding  to 
conjugate  complex  roots  X3  and  X4,  radians 
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9 


_  <^n 

np "  aF 


&0>al>a2>a3>ak>&ti)  coefficients  appearing  in  numerator  terms  of  ampli¬ 
tude  coefficients  for  0 

coefficients  appearing  in  numerator  terms  of  ampli- 
r  tude  coefficients  for  \|r 

c0tcl>c2>c3>ck  coefficients  appearing  in  numerator  terms  of  ampli¬ 

tude  coefficients  for  p 


Subscripts: 

0  initial  value 

0  transformed  variable 


ANALYSIS 


The  linear  equations  of  motion,  referred  to  the  axis  system  shown 
in  figure  1  and  representing  the  lateral  motion  of  an  airplane  are 


10 


NACA  TN  2129 


2Mt)Kx2l>b20  "  JClpM  +  2^KXZDb2llf  "  "  CZpP  “  C2c 

2nbKxzDb2^  -  ic^DbSZf  +  2^Kz2Dt2i  -  -  CnpP  - 

-  CL0  +  2^bV  -  CL  tan  7*  .  -  CYpP  + 


2MbDbP 


The  terms  C7  ,  Cn  ,  and  Cy„  are  forcing  functions  which  represent 
c  c  c 

disturbances  imposed  upon  the  state  of  motion  of  the  airplane  by  control 
movement  or  atmospheric  turbulence.  These  terms,  in  general,  are  arbi¬ 
trary  functions  of  time,  but  for  the  purpose  of  this  analysis,  they  are 
considered  to  be  constants  applied  at  zero  time.  After  a  solution  has 
been  obtained  in  terms  of  constant  forcing  quantities  this  solution  can 
be  used  to  obtain  a  new  solution  for  an  arbitrary  forcing  function  by 
Duhamel's  integral  as  explained  in  references  4  and  5* 


Transformation  of  Equations 

'When  the  Laplace  transform  is  applied  (reference  6,  p.  8),  the 
transformed  equations  become  after  multiplying  through  by  0 


(2HbKx2ff3  ‘  !%a2H  +  (2^KXZa3  “  lflra2)^o  +  ('Czpa)PcJ  =  rl 


rl  =  (2MbKxV  -  |%a)0o  +  (2ubKxza2  -  §CZrcr)*0 
(2|ibKx2cx)(l>b0)o  +  (2^^)  (hb^Q  +  CZ(; 


\  (2a) 
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(WjE®3  -  |CV2)^  +  (^Z8®3  -  !<V2K  +  (-CnjOjPa  -  C 

r2  =  (2^KXZa2  "  |Cnp^0o  +  (2^z2°2  "  gCnra)  ^0  +  /  (2b) 

(2%KXZa)  (l>b0)o  +  (2MbKz2«y)(Db+)0  +  Cnc 


(-|-cYpCT2  -  C '■Lo)$a  +  (2^a2  -  ^a2  -  CL  tan  7cr)ta  + 
(2^a2  -  CYpo-)3a  =  r3 

r3  =  (-FCYpa)^0  +  (2^a  -  ^YPff)t0  +  (2%°K  +  CYc 
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and  the  constants  are  given  by 


a0 


%  V  c*r 


=  0nkCnR  C.  KZ2  +  ldc 


CZ„  KX2 


nR  ^XZ  ^n. 


2  ‘P 


CY  Cy  0 
*|3  P 


CYp  0  CY; 


Cnp  KXZ 
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% 

04 

X 

X 

Kxz 

D  -  2ih,CL  Cjjp 

Kxz 

Kz2 

0 

l 

tan  ; 

cl 

% 

C2r 

K  -  -$°i  % 

% 

cnr 

0 

1 

tan  7 

%  %  V 

l  cn  cn  Cnr  Hr  Hb  *P  CZp 

Cnp  Cn 

%  %°lr 


The  expression  for  is 

■ir  _  bO<*5  +  bicA  +  b2cr3  +  b +  b^cr  +  b^ 

Yg-  —  -  ” 

ct2A 

where  the  constants  are  given  by 


bo  =  *b  P^b3 


Ky  K, 


Kxz  Kz2 


/  KX2  KXZ 

+oRb2Kxz  kz2 


2nb2  ClP  ^  +  (DfetdsMb3  Kx  KXZ 

Cnp  KZ  l  %z  Kz2 
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H  =  iUj  cL 


^np  Cnp 


+  d>b0)o  A^cL|CZp  Kx2 


ICnp  K 


XZ 


0 


/ 

CZR  KXZ 

%  %  Czc 

|_2MhcL 

P  ^ 

Cnp  Kz2 

)•> 

Cnp  cnp  CDc 
CYp  CYp  CYc 

b5  =  -CL 


cz  C7 
4p  lc 

°np  Cnc 


The  expression  for  $n  is 


Pa  = 


c0cA  +  c-^a3  4-  c^cr2  4-  c^cr  4-  cj^ 


crA 


where  the  constants  are  given  by 


c0  =  Pol8^3 


K> 


XZ 
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C1  =  00  I^cl 


KX  Kxz 


Kxz  Kz' 


+  (Db0)n  2^2Cy 


Kxd  KXZ 

Kxz  Kz2 


Kxz  Kz2 


I  Kx2  Kxz  \  ( 

^b2tan  7  CL  +  (Db^)0  2Mb2  (cYr  -  ^Mb) 


KX2  KXZ 


KXZ  Kz' 


a  L  2  °lr  KX2  n  p  %  KXZ|\  ,  2  KX2 

P0  2Mb2  -  2Mb2  P  J  +  H^Y, 


Cnr  Kxz 


%  KZ2 


Kxz  Kz2 


=  0o(^bCL 


ci  Kx2  civ  kyz\  I  Clp  Clr 

/  -PbCL  P  UcW^bKx,  Cnp  Cnr 

°nr  Kxz  Cn  Kz*  /  l 

/  \  0  %  % 


‘r  iVXZ 


KX  CZr,  KXZ 


C,  K-. 


1  /  ^  Z  1VX 

2ha2  KXZ  cnD  Kz2  +  to  Mb  tan  7Cl  r 

\  C%  Kxz 

0  2CL  1  /  \ 


lCZ-D  Kxz 


Cnp  KZ‘ 


Mb  tan  7Cl  *  .  ]  +  (E>b^)o  i^b  KZ2  cn-n  Cnr  + 

n  v  d  c  JP 


KXZ  C*t)  °lT 


0  cYp  cYr 


2^  K 


KZ2  +P0/kCZp  'HU 


1  -2  tan  7Cl  0 


2nr 


0  KX^  Ke 


Czc  %  cZr 


Cnc  cnp  cnr  -  ^b  0  KXZ  Kz2  -  ^ 


Cy  Cy  Cy 
c  p  r 


Cy  Cy  Cy 
±c  xp  xr 


2  Kx2  Clc 
Kxz  Cnc 
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All  the  determinants  given  in  this  paper  are  expanded  in  the 
appendix. 

In  order  to  obtain  the  actual  variables  0,  \|r,  and  |3  from  the 

transformed  variables  an  inverse  LaPlace  transformation  must  be  applied 
t0  $a)  and  3a.  The  expressions  for  <fa,  and  j3a  are  of  the 

form  p a/%  where  pCT  and  qa  are  polynomials ,  the  degree  of 
being  higher  than  that  of  pa.  Reference  6,  page  4 5,  indicates  that  the 
inverse  transform  of  a  function  of  this  type  is  (in  terms  of  the  variables 
used  herein) 

I^-l/PgN  _  Pg(  ^n)  e  *-nsb 

W  ~  Si  V(xn)  6 

This  equation  assumes  all  the  roots  \n  of  q0  =  0  to  be  distinct.  All 

roots  of  q0  =  0  are  distinct  for  however,  for  (fa  and 

q<j  =  0  has  double  zero  roots.  (See  equations  (3),  (5),  and  (6).)  The 
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terms  in  the  equations  for  0  and  >Jr  resulting  from  the  two  zero  roots 
of  =  0  are  given  according  to  reference  6 ,  page  ^9>  hy 


— (0)  +  fi(0)sb 

da 

where 


ft  = 

% 


The  analysis  takes  three  forms  depending  upon  the  character  of  the 
nonzero  roots  of  q^  =  0  which  are  the  same  as  the  roots  of  A  =  0. 
Four  real  roots,  two  real  roots  plus  a  pair  of  conjugate  complex  roots, 
or  two  pairs  of  conjugate  complex  roots  may  exist. 


Four  Real  Roots 

The  inverse  Laplace  transform  of  is 


0  =  Axe 


xlsb 


+  A^e 


X2sb 


+  A^e 


XoS 


3bh 


A^e 


x4sb 


A58b  +  Ag 


(7) 


where  the  constants  are 


Ai  = 


a.r\Xn5  +  St  Xt  ^  3  +  StX-,2  +  a^X^  +  a, 


‘■err 


ii 


2  1 


X_i 


6aXx5  +  5BX-L11  +  hCX-^  +  3DX12  +  2EXX 


a0X25  +  a^X^  +  a2X2^  +  +  al*X2  +  &5 

h-2  =  ■  -  "  ■  - 1  -  - - ' 

6AX25  +  52X2^  +  1*CX23  +  3DXg2  +  2EX2 


A3  = 


a0x3^  +  alx3^  +  a2x3^  +  a3x32  +  a4x3  +  a5 
6a  x35  +  5bx3a  +  lie  x33  +  3dx32  +  2EX3 
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A  _  aQ^  +  +  a-2^3  +  a^X^2  +  a^X^  +  a. 

6A\ip  +  [jBXi^  +  ItCXjj3  +  3DX42  +  2EX4 


A=  =  S 

5  E 


H  -  |(h  -  »s|) 


The  inverse  Laplace  transform  of  is 


♦  *  B^1Sb  *  V8*"  ♦  B^3*»  +  ^  +  ^  + 


where  the  constants  are 


»!  -  ^ 1  b^ig  *  b^i  j  h 

6AX +  5BX1li’  +  teX-j^  +  3BXX2  +  2EX-L 


^0^2^  +  ^1^2^  +  +  b^Xg2  +  bj^Xg  +  b, 

6a^5  +  58X2^  +  1*CX23  +  3DX22  +  2EI2 


=  V"3^  +  blx34  +  b2X33  +  b3^32  +  +  b5 

6aX35  +  5BX2^  +  1»CX33  +  3DXo2  +  2EXo 


^  +  bl^J+  +  b2X43  +  b^X^2  +  bi^  +  h^ 
6axi,.5  +  5BX^  +  iicx^3  +  30x^2  +  2EX4 
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B6  -  Kk  -  t>5§) 

The  inverse  Laplace  transform  of  (3a  is 

P  =  CieXlS^  +  Cge^b  +  C3e^3sb  +  c^k°b  +  c5 

where  the  constants  are 


Cl  = 


cQXl^  +  cix-jl4'  +  C2X1^  +  C3X12  +  C4X1 
6AXx5  +  5BX.11+  +  kCX^  +  3D^X2  +  2E^x 


C2  = 


c0^2^  ci^2^  c 2^2^  ^  c 

6kXp  +  5B\2*  +  bCX^  +  3DX22  +  2E^2 


=  cox3^  +  cix3^  +  C2X33  +  C3X32  +  c4X3 
GkX-p  +  5BX.3^  +  i4CX33  +  3DX32  +  2EX.3 


CJ+ 


=  cQX^  +  cx^  +  c2X43  +  c3X42  + 

Gk\\p  +  5BX-4^  +  itCX^  +  3Dx.^2  +  2EX4 


Cr  = 


ch- 


(9) 
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The  quantity  can  he  obtained  from  equation  (7)  by  differentiation 


Db0  =  Ai'e^  +  A2'e^Sb  +  A,'eX3*b  +  Ak'e^  + 


where  the  constants  are 


—  AqA}. 


Ag!  -  X2A2 


A3'  _ 

A^'  =  X^A^ 

V  =  a5 

The  quantity  can  be  obtained  from  equation  (8)  by  differentiation 

as 

Dbt  =  B1,eXl8b  +  B2'eX2Sb  +  B3'eX3&b  +  B4'eX^Sb  +  B5'  (11] 

where  the  constants  are 


— 


=  X4B4 
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The  quantity  D^P  can  be  obtained  from  equation  (9)  by  differentiation 
as 

D„0  =  +  C2'eX2Sb  +  C3'eX3Sb  +  Cl(.eH=b  (l2) 

where  the  constants  are 

cl!  =  Xlcl 
C2 T  =  ^2^2 
c3’  =  X,3c3  . 

C4  —  X4C4 


Collecting  the  equations  of  motion  (equations  (7)  to  (12))  for  the  case 
of  four  real  roots  gives 

i  =  A^1^  +  A2eX28b  *  A3eX3S»  +  A,,.'*’6  +  A,^  +  A^ 

*  =  BieXlSb  +  B2eX2S1>  +  B3eX3sb  +  B^1*^  +  +  B6 

B  -■C1«Xl’b  +  C2eX2S6  +  C3eX3sb  +  +  c? 


\  (13) 

=  A^e^  +  A2-eX2Sb  +  A^e^  +  A^'e^  +  A y 

Db*  =  B1'eXlSb  +  B2'eX2Sb  +  B3'eX3Sb  +  B^'e*^  +  B 5' 

DbP  =  C1'eXlSb  +  C2'eX2Sb  +  C3'eX3Sb  +  Cj/e^h  ^ 

Two  Real  Roots  and  a  Pair  of  Conjugate  Complex  Roots 


If  a  pair  of  conjugate  complex  roots  X.3  and  exists,  the 

coefficients  of  the  terms  of  0,  t,  and  3  (equation  (13))  corresponding 
to  A.  3  and  XI4.  are  conjugate  complex.  The  complex  number  X3k  can  be 
written 

\3k  =  Rk  +  Iki 


NACA  TN  2129 


25 


Thus , 

*3  =  Rl  +  Iii 
=  ^2  +  IqI 

i^i 

l3k  =  R|,  +  I|,i 
X35  =  Bj  +  I,i 

The  coefficients  of  the  terms  of  0  resulting  from  the  complex  roots 
are 
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The  terms  of  0  corresponding  to  the  conjugate  complex  roots  are 
A  X33b  +  A^e^sh  _  KAeRlSb  costly  +  u>A) 


where 


and 


K, 


=  2  \Jba2  +  IA2 


.  -1  lA 

“A  =  tan  / 


Similarly,  these  terms  for  t  are 

B^e^8*3  +  =  KBeRlS^  cos(lisfc,  +  cob) 

where 


and 


KB  =  2\Jrb2  +  IB2 


o>B  =  tan"1  ^2 
rB 


and  for  (3  are 


Coe^3Sb  +  c^e^8*5  =  KceRl8b  costl-^  +  o^-.) 


where 


and 


K, 


Rc2  +  Ic2 


-1  Tc 

=  tan 
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Two  Pairs  of  Conjugate  Complex  Roots 

If  two  pairs  of  conjugate  complex  roots,  Xi,  and  X3, 

exist,  another  cosine  term  is  introduced  into  equations  (l4)  in  place 
of  the  exponentials  so  that 


0  =  KA'eRl  cos(l1’sb  +  mA') 
t  =  KB'eRl  Sb  cos^I^'Sh  +  u>g'^ 


+  KAeRl&b  cos  ^I]_sb  +  cjda^  +  +  Ag 


~\ 


+  Kge^313  cosp1sb  +  cogj  +  B5sb  +  B6 


P  =  KQ'eRl  Bl°  cos^Ib'sb  +  +  KQeR-*-&t-’  cos^Ibsb  +  +  C5 


°b0  =  kA '  \Jft i'2  +  ^l'2  e  1  ^  co^Ib ' 


.  +  -1 TA 

sb  +  ct>A '  +  tan  -  + 

R1V 


+  Ib2  eR^Sl3  cos  ^bsb  +  cda  +  tan-1  +  A^ 

Dbt  =  V^V2  +  ^  %  c0B(^l'3-b  +  V  +  tan_1  ^7^  + 

KB  eRl8^  cos  ^Sh  +  0>B  +  tan"1  h  'j  +  B 5 

DhP  =  +  Ib'2  e  13  cos^Ih'Sh  +  o^-.'  +  tan"-1-  + 


K, 


\/Rl2  +  *1 


Q  Rl  S-U  /  1  Ii 

^  e  ^  u  cos  !]_%  +  cDq  +  tan"  — 


>(15) 
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4 

DISCUSSION 


The  lateral  motions  of  aircraft  have  been  obtained  by  means  of  the 
Laplace  transform.  This  analysis  resulted  in  equations  from  which  the 
free  lateral  motion  of  an  aircraft  can  be  calculated  for  any  initial 
condition,  or  the  forced  motion  can  be  calculated  for  any  constant  lateral 
force  or  moment  applied  at  zero  time.  In  general,  the  lateral  forces  and 
moments  applied  to  the  airplane  by  control  movement  or  atmospheric  tur¬ 
bulence  are  not  constant  but  are  arbitrary  functions  of  time.  After  a 
solution  has  been  obtained  in  terms  of  constant  disturbing  forces  and 
moments,  however,  the  solution  for  the  arbitrary  forces  and  moments  can 
be  obtained  by  Duhamel's  integral  as  explained  in  references  4  and  5. 

The  nature  of  the  motion  indicated  by  equations  (13)  to  (15)  depends 
upon  the  form  of  the  roots  of  the  polynomial 

A  =  0  ' 

which  is  commonly  referred  to  as  the  stability  quartic.  The  roots  of  the 
quartic  can  take  three  forms  -  first,  all  four  roots  real)  second,  two 
real  roots  and  a  pair  of  conjugate  complex  roots)  and  third,  two  pairs  of 
conjugate  complex  roots.  In  the  case  of  the  lateral  motions  of  airplanes 
the  first  form  almost  never  occurs)  the  second  form  is  very  common)  and 
the  third  form  occurs  under  rather  rare  conditions.  The  actual  motions 
indicated  by  equations  (13)  to  (15)  can  be  seen  to  be  composed,  in-  general, 
of  the  sums  of  terms  which  are  the  amplitude  coefficients  (the  A's,  B's, 

C  s,  and  K  s  of  equations  (13)  to  (15))  modulated  by  exponential  and  cosine 
factors. 


All  the  classical  stability  concepts  can  be  obtained  from  equations 
(13)  to  (15).  Because  stability  is  concerned  only  with  the  free  motion 
(motion  due  to  initial  conditions)  the  forcing  or  disturbing  quantities 
^lc>  Cnc>  an<l  CYC  can  be  set  equal  to  zero  so  that  the  amplitude  coef¬ 
ficients  Aij,  Be;,  and  C<j  vanish.  The  variation  of  the  amplitude  of  the 
motion  with  time,  which  determines  the  stability,  is  now  dependent  entirely 

upon  the  damping  coefficients  e^1&b,  e^2&b,  e^8!  e^4sb  eRlsb 

R  1  «  7  7 

and  e  1  The  motion  diminishes  with  time  (stable)  if  X^,  X2,  X3, 

X4,  R]_,  and  R]_'  are  all  negative.  Thus,  these  criteria  for  stability 
are  that  all  real  roots  of  A  =0  and  the  real  parts  of  all  complex 
roots  be  negative.  These  criteria  have  been  expressed  in  reference  7  in 
terms  of  the  signs  of  the  coefficients  of  the  quartic  A  =  0  and  the 
sign  of  Routh's  discriminant  which  is  written  as 

R  =  BCD  -  AD2  -  EB2 


(16) 
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These  criteria  in  the  present  case  can  he  expressed  as  follows:  The 
necessary  and  sufficient  conditions  for  the  real  roots  and  the  real  parts 
of  the  complex  roots  to  he  negative  are  that  every  coefficient  of  the 
quartic  and  R  should  he  positive. 


If  the  motions  contain  oscillations,  the  periods  of  the  oscillations 
in  seconds  are  from  equations  (ll)  and  (15) 


P 


2  rib 


X1V 


P' 


2ffb 

Il’V 


(17) 


and  the  times  to  damp  to  half-amplitude  in  seconds  are 


Tl/2 


’1/2 


h  loge2 
r-lV 


,h  loge2 
r-l'v 


~\ 


J 


and  the  cycles  to  damp  to  half -amplitude  are 

„  ^l/2  ^1  ^°®e^ 

»i/2  ■  -i — 


t  _  ^l/2  _  Ij  * 

l/2  "  "pT  2«Rl' 


(18) 


(19) 


APPLICATION 


The  equations  for  the  motion  of  an  aircraft  resulting  from  the 
analysis  were  used  to  calculate  illustrative  examples  of  certain  free 
and  forced  motions  of  an  experimental  swept-wing  airplane,  a  three-view 
drawing  of  which  is  shown  in  figure  2.  The  calculations  were  made  hy 
use  of  the  Bell  Telephone  Laboratories  X-667H  relay  computer  available 
at  the  Langley  Laboratory.  The  calculations  were  based  upon  stability 
derivatives  measured  on  a  model  of  the  experimental  airplane  in  the 
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Langley  stability  tunnel  and  presented  in  reference  8.  Motions  were  cal¬ 
culated  for  true  airspeeds  of  140  and  200  miles  per  hour  under  standard 
conditions.  The  stability  derivatives ^  other  related  aerodynamic  quan¬ 
tities,  and  the  mass  characteristics  of  the  experimental  airplane  as  used 
in  the  calculations  are  shown  in  table  I.  The  coefficients  of  the  sta- 
b il i ty  quartic  (equation  (4))  and  value  of  Routh's  discriminant  (equa¬ 
tion  (16) )  were  calculated  as  * 


CL 

True 

airspeed 

(mph) 

A 

B 

c 

D 

E 

R 

0.693 

140 

26.19791 

10.18804 

3.021074 

0.6312249 

0.00235618 

5.8 

•  340 

200 

26.20030 

9.818377 

cS 

o- 

S' 

LP\ 

OJ 

.4623735 

.00014875 

8.7 

The  positive  signs  of  all  these  quantities  indicate  complete  stability 
of  the  lateral  motion  for  both  airspeeds.  The  coefficients  of  the  quartic 
are  such  as  to  give  two  real  roots  and  a  pair  of  conjugate  complex  roots 
which  are 


True 

cL 

airspeed 

(mph) 

X1 

X2 

E±  ±  il± 

0.693 

140 

-0.2802853 

-0.003603100 

-0.05249952  ±  0.2859079i 

•  340 

200 

-.2649690 

-.0003222716 

-.05472583  ±  0.2519754i 

Discussions  of  methods  of  obtaining  the  roots  of  the  quartic  can  be  found 
in  references  9  to  12, 


The  first  or  second  powers  of  o*  multiplied  by  A  in  the  denomi¬ 
nators  of  equations  (3)*  (5)>  and  (6)  introduce  one  or  two  zero  roots, 
respectively,  in  addition  to  the  roots  given  in  the  preceding  table. 

These  zero  roots  lead  to  the  terms  containing  the  amplitude  coefficients 
Ag,  Bg,  and  C5  of  equations  (l4).  For  the  experimental  air¬ 

plane,  the  motion  can  be  thought  of  as  composed  of  three  modes  -  the 
oscillatory  mode  resulting  from  the  pair  of  conjugate  complex  roots,  the 
rolling-subsidence  mode  resulting  from  the  large  negative  real  root  Xi 
and  the  spiral  mode  resulting  from  the  small  negative  real  root  and  the 
zero  roots.  Stability  of  the  free  spiral  motion  is  indicated  by  the 
negative  sign  of  the  small  real  root  X2. 
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The  period,  time  to  damp  to  half -amplitude,  and  cycles  to  damp  to 
half -amplitude  of  the  oscillatory  motion  were  calculated  by  use  of  equa¬ 
tions,  (17)  to  (19)  from  the  imaginary  and  real  parts  of  the  conjugate 
complex  roots  as 


CL 

True  airspeed 
(mph) 

P 

(sec) 

Tl/2 
( sec) 

Nl/2 

0.693 

140 

3.60 

2.16 

0.60 

.340 

200 

2.86 

1.45 

.51 

The  motions  calculated  fell  into  two  categories  which  may  be  termed 
free  and  forced  motions. 


Free  Motions 

Free  motions  are  those  which  exist  following  an  initial  condition 
and  in  the  absence  of  any  forcing  function.  The  five  possible  initial 
conditions  are  0q,  \|Tq,  |3q,  rQ,  and  pQ.  Every  free  motion  that  the 

airplane  is  capable  of  executing  can  be  obtained  by  superposition  of  the 
motions  following  these  initial  conditions  taken  separately.  Figures  3 
to  6  show  the  calculated  free  motion  following  the  initial  conditions  0q, 
(3q ,  tq,  and  pq  for  the  experimental  airplane  in  level  flight  according 
to  equations  (l4).  No  airplane  response  to  occurs  when  the  angle  of 

climb  is  zero.  Table  II  gives  the  values  of  the  amplitude  coefficients 
(see  equations  (l4))  corresponding  to  the  motions  of  figures  3  to  6.  These 
figures  show  the  total  airplane  motion  and  show  the  separate  contributions 
to  the  motion  by  the  rolling-subsidence  and  spiral  modes  when  the  motion 
resulting  from  these  modes  is  appreciable.  Figures  3  to  6  indicate  that 
in  the  case  of  the  experimental  airplane  the  initial  condition  (3q  pre¬ 
dominately  excites  the  oscillatory  mode  of  motion,  tq  excites  both  the 
oscillatory  and  spiral  modes,  and  and  pq  predominately  excite  the 

spiral  mode.  The  rolling- subsidence  mode  appears  for  a  very  short  period 
of  time  in  the  initial  phases  of  any  motion  involving  appreciable  rolling 
velocity.  The  principal  effects  upon  the  motions  of  figures  3  to  6  of 
increasing  the  airspeed  from  140  to  200  miles  per  hour  is  to  reduce  the 
period  as  well  as  the  time  and  cycles  to  damp  to  half- amplitude. 


Forced  Motions 

Forced  motions  are  those  which  exist  during  the  action  of  forcing 
functions  upon  the  airplane.  Any  forced  motion  of  an  airplane  can  be 
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built  up  by  proper  superposition  (Duhamel's  integral,  reference  4)  of  the 
motions  following  the  application  at  zero  time  of  constant  values  of  the 
forcing  functions  C 7  ,  Cn  ,  and  Cy  .  Figures  7  and  8  show  the  cal- 

°  C  C  # 

culated  response  according  to  equations  (l4)  of  the  experimental  airplahe 

to  the  constant  value  0.02  for  Ci  and  Cn  applied  at  zero  time.  The 

C  X1c 

response  to  a  value  of  0.02  for  C yc  was  also  calculated  but  during  a 
time  period  of  8  seconds  was  negligible  compared  with  responses  resulting 
from  Cic  and  Cnc.  For  the  experimental  airplane,  the  value  C 1  =  0.02 

corresponds  to  a  total  aileron  deflection  of  21.0°,  the  value  of  Cnc  =  O.i 
corresponds  to  a  rudder  deflection  of  13.7°,  and  the  value  Cy  =0.02 
corresponds  to  a  rudder  deflection  of  7-5°.  The  response  to  C 1  is  pre¬ 
dominately  in  the  spiral  mode  of  motion;  whereas  the  response  to  Cn  is 

nG 

predominately  in  the  spiral  and  oscillatory  modes. 

A  large  number  of  forced  motions  calculated  for  the  experimental  air¬ 
plane  corresponding  to  various  flight  conditions  are  presented  in  refer¬ 
ence  8.  These  motions  were  built  up  by  superposition  of  motions  such  as 
those  of  figures  7  and  8  following  the  application  of  the  constant  forcing 

functions  Cj  and  Cn  .  A  large  number  of  comparisons  are  made  in 
c  c 

reference  8  between  calculations  and  flight  tests  for  a  large  variety  of 
flight  conditions.  The  agreement  between  calculated  and  flight  motions 
is  good  and  indicates  the  practicability  of  analyzing  the  dynamic  lateral 
qualities  of  aircraft  by  use  of  the  theory  of  lateral  dynamics  such 
as  that  herein  developed  if  experimentally  determined  values  of  the  aero¬ 
dynamic  and  mass  parameters  of  the  airplane  are  available.  Reference  8 
also  indicates  rather  insignificant  effects  upon  the  calculated  motions 
that  result  from  a  consideration  of  slight  nonlinearities  which  occur  in 
certain  aerodynamic  parameters  of  the  experimental  airplane. 


CONCLUDING  REMARKS 


The  lateral  motions  of  aircraft  were  determined  by  means  of  the 
Laplace  transform  which  gave  solutions  expressed  in  terms  of  elementary 
functions  for  the  free  and  forced  motions.  These  equations  permit  the 
calculation  of  the  free  motion  of  an  aircraft  following  any  initial  con¬ 
dition  or  the  forced  motion  following  the  application  of  constant  external 
forces  and  moments.  These  forced  motions  can  be  used  to  obtain  the 
response  to  any  arbitrary  forcing  function  by  means  of  Duhamel's  integral. 
All  the  classical  stability  concepts  can  be  deduced  from  these  same 
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solutions  largely  by  inspection.  These  equations  for  the  lateral  motion 
were  applied  to  the  calculation  of  the  lateral  stability  of  a  specific 
airplane  and  to  the  calculation  of  certain  of  its  free  and  forced  motions. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  April  3>  1950 
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APPENDIX 

EXPANSION  OF  THE  COEFFICIENT  DETERMINANTS 
The  stability  quartic  coefficients  are 
A  =  8M.-b^Kx2Kz2  -  Kj^2) 

B  .  -  KX'%2)  +  Kxz(C,r  +  C„p)  -  Kx2Cnr  -  KZ2C,1 

c  -  V(kx2%  -  *xz%)  +  %%2(c%cY|3  -  cnpcIr)  +  - 

%%)  +  ^KXZ  (%%  -  CV*P  +  %%  -  %%)  + 

(C 2pCnr  "  c*rCnp) 


E 


=  2^CL|tan  7^ CZp  -  KX2cnp)  +  (KXZCnp  “  KZ2c2p)  +  (CZpcnp  " 
cZpCnp),+  ^(c2rCnpCYp  +  Cl^Cn^Cjr  +  CzpCnrCYp  -  C ZpCnpCyr  - 

C2pCnrCYp  -  c2rcnpcYpj 

*  |cL(Ian  r(c,pCnp  -  C,BCnp)  +  CI@C„r  -  C!rcJ 
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The  coefficients  appearing  in  the  numerator  of  amplitude  coeffi¬ 
cients  for  0  are 

ao  =  0oJ^3(kX2kZ2  -  kXZ2)J 

ai  =  0o[2p-b2Kxz(c7r  +  Cnp)  +  H2cYp(KXZ2  '  ^Z2)  '  2^b2(KX2cnr  + 

KZ2c2p|  +  (^)o  8Mh3(KX2KZ2  "  KXZ2)] 
a2  =  00  ^Kx2(CnrCYp  -  CnpCYr)  +  ^2(KX2Cnp  -  KXZCzp)  +  M.bKZ2  (%%" 
clpcYp)  +  RjKXZ  (CnpCYp  -  CnpCYp  +  CzpCYr-  C2rCyp^  +  |nb(C2pCnr  - 
C2rCnp)|  +  (Db0)o  ^Mb2CYp(KxZ2  -  Kx^Z2)  +  2Pb2(KxzClr  -  KX2Cnr)j  + 
(Db+)o  2Mb2^KZ2Czr  -  KxzCnr)l  +  PO  Vb2(KZ2Czp  -  KxzCnp)  + 

^Mb2  (kZ2cZc  '  KXZCnc) 


-MM 
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a4  =  0o[|  tan  7CL  (czpCn3  -  CZpCnp)  +  (^)o  [2%  tan  ?CL  (Kxz%- 

KX2cnp)j  +  +o  \  tan  7  CL(c2rCnp  -  C^C^j]  + 

(Db^O^t,  tan  7CL  (kz2C2|3  -  K^C^J  +  + 

lClc(cnrCYp  -  CnpCYr)  +  \  Cnc(cipCYr  -  CzrCYp)  + 

|  CYc(czrCnp  -  CzpCnr)J 
a5  =  tan  7CL  (C2pCnc  -  C2c%) 

The  coefficients  appearing  in  the  numerator  of  amplitude  coefficients 
for  i  are 

b0  =  to  8^b3  (kx2kz2  -  Kxz2)] 

bX  =  Mr0  2nb2Kxz(cnp  +  Czr)  -  ^2CYp(Kx2Kz2  -  *XZ2)  -  2^2(KZ2CZp  + 
KX2Cnr)]  +  (^t)0[8^3(KX2Kz2  "  KXZ2)] 
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*2  =  M)o^b2(KX2cnp  “  KXZ%)]  +  ^o[^KX2(CnrcYp-  CnpCYr)  + 

^Kz2(c2pCYp  -  CZpCyJ  +  ftKxz  (cnpCyp  -  CnpCyp  +  CZpCyr  -  CzrCyp)  + 

^2(KX2cnp  '  KXZcZp)  +  \  **b(clpcnr"  C*rCnpj]  + 

(Db^)o[2^2(KXZCnp  -  Kz2Czp)  -  ^b2Cyp(KX2KZ2  -  + 

Po[K2(KX2Cnp  -  K^C^I  +  H*2(Kx2cnc  '  KXZCIc) 

b3  =  0oj2^bcL(KX2cnP  -  KXZCZPJ  +  (^f^P  "  %%)  + 

MbKxz(ClpOlp  '  %%)]  +  *o[^(CZp%  “  %%)  ‘  " 

KXZcnp)  +  ^(cZrCnpCYp  +  cZpcnpcYr  +  cZpcn1.Cyp  -  cZpcnpCYr 
Czp<V*P  -  CzrCnpCYpj]  +  (l>b*)o[w*Z8 (czpCyp  -  Clp(*p)  + 

^KXz(cnpCYp  “  %%)]  +  Po[^b(CZpCHp  “  %%)] 

^(Czc%  -  CZpCnc)  +  2^E(clcCjrp  -  Cl(jCYc)  + 

2^KX2(GnpCYc  "  CncCYp) 


ko 
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H  =  ^oj^L^ZpCnp  -  ClpCnp|  +  (Db0)o  j2^CL^Kx2Cnp  -  KxzCZpjJ  + 
t0|cL(czpCnr  -  CzrCnpj]  +  (Dbt^^bCL^xzCnp  -  KZ2Czpjj  + 
|czc(cnpCYp  -  CnpCYpj  +  Jpnjci^  ~  %%)  + 

!%(%%  -  %%) 
b5  =  cL(czccnp  -  c2pcnc) 


The  coefficients  appearing  in  the  numerator  of  amplitude  coeffi¬ 
cients  for  3  are 


co  =  3o[q^3(kx2Kz2  -  Kxz2] 

Cl  =  0o[^132cl(KX2kZ2  -  KXZ2)] 

yQ  L*2  tan  rCL  (kX2Kz2  - 


+  (l>b0)o  2^b2cYp(KX2KZ2 
Kxz2)j 


Kxz2)]  +  P0[^2(KXZCZr  -  Kx2Cnr)  +  2^2(KXzCnp  -  KZ2Czpjj  + 

Vb2CYc(Kx2Kz2  -  Kje2) 
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c2  =  0o[^bcL(KxzCzr  -  Kx2cnr)  +  Mb'CL^C^  -  + 

(ciipCYr  -  CnrCYp)  +  2[i^X2^KZ2CL  -  Cn^  + 
2%2%z(cZp  ‘  ^XzCl)  +  HKXz(ClrcYp  -  CZpCYrj]  + 

♦0[*  tan  7CL  (KXZcZr  -  C^Kx2)  +  Hb  tan  7CL  -  Kz2C2jJ  + 

(Db^)o[i^KXz(CnpCYr  -  C^CyJ  +  |^Kz2^2rCYp  -  C2pCYr)  + 

H2  tan  7CL  (kx2Kz2  -  K^)  +  2,^Kz2CZp  -  K^C^J  + 

P°[iU'b(CzpCllr  "  CzrCnpj]  4  4^2(kxzCzc  -  KX2Cnc)  +  ^KZ2^cCYp  - 

C2pCYc)  4  ^KX2(CncCYr  -  CnrCYc)  +  ^bKxz^CnpCYc  -  + 

^KXZ^c2r°Yc  "  CZcCYr^ 
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c3  =  -  C^)]  +  M^L  tan  7(KX2Cnp  -  KxzPlp)  + 

^bCL^xzCZr  -  KX2Cnrjj  +  *o[|  tan  7CL  (C*pCnr  “  C*rCnP)]  + 

(^)0[^CL(Kz2CIr  -  KxzC^)  +  Ht,  tan  7CL  + 

^(%cnc  -  CZcC Dp)  +  2^cL  (kZ2cZc  '  KXZcnc)  + 

2Mb  tan  7CL  (Kx2Cnc  -  K xzcZc)  +  j(CZcCnpCYr  +  ClpCnrCYc  + 
cZrCncCYp  "  CZr%CYc  "  %CncCYr  '  C*cC*r%) 

=  ^lrCnc  -  ClcCnr)  +  2  7°L  KS  *  %Cnc) 
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TABLE  I 

AERODYNAMIC  AND  MASS  CHARACTERISTICS  OF  AN 
EXPERIMENTAL  SWEPT -WING  AIRPLANE 


l4o 

200 

w 

8700 

8700 

s 

250 

250 

m 

270.2 

270.2 

P 

0.00238 

0.00238 

7 

0 

0 

1) 

33-6 

33.6 

CL 

0.693 

0.340 

a 

10.6 

4.8 

n 

11.05 

5.25 

13.51 

13.51 

v/b 

6.111 

8.730 

KX2 

0.02329 

0.02219 

OJ 

IS] 

0.05932 

0.06042 

Kxz 

0.007316 

0.003544 

% 

-0.0659 

-0.0275 

•  Cnp 

0.100 

0.0975 

% 

-0.739 

-0.722 

% 

-0.325 

-0.31 

% 

-0.1 

-0.06 

Cy 

1P 

0.44 

0.3 

c2 

Lr 

0.12 

0.07 

Cn 

nr 

-0.280 

-0.280 

Cy 

xr 

O.36 

0. 40 

Hliles-per-hour  units. 
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TABLE  II 

AMPLITUDE  COEFFICIENTS  OF  FBEE  AND  FORCED  MOTIONS  OF  EXPERIMENTAL  SWEPT-WING  AIRPLANE 


(a)  V  =  140  miles  per  hour. 


Initial  conditions 
and  forcing 
functions 

Variables 

Rollinu-BtiDsidence 

mode 

Oscillatory  mode 

Spiral  mode 

0 

A1 

0.0 4073926 

% 

0.05404332 

a2 

0.4374647 

a5 

0 

A6 

0 

+ 

Bl 

-.00222650 

kb 

.04009448 

b2 

-3.038911 

b5 

0 

b6 

3.029296 

*0 

P 

Cl 

-.00131258 

Kc 

.04330260 

C2 

.01392006 

c5 

0 

- 

P 

Al' 

-.06978088 

KA\/Rl2  +  Il2 

.09600416 

V 

-.00963249 

V 

0 

■r 

V 

.00381366 

Kb^I2*1!2 

.07122481 

b2’ 

.06691349 

b5* 

0 

- 

0 

Al 

-.01808626 

Ka 

.2450096 

a2 

-.02458282 

a5 

0 

a6 

0 

If 

B1 

.00973284 

Kb 

.18177064 

B2 

.17076788 

b5 

0 

B6 

0 

Po 

P 

Cl 

.00573756 

Kc 

.19631484 

C2 

.00078222 

c5 

0 

-- 

P 

Ai* 

.30503^2 

>WBia  +ii2 

.43524085 

A2 

.00054129 

a5' 

0 

- 

r 

Bi' 

-.01667089 

KB^12  + 1]_2 

.3229020 

V 

-.00376012 

b5' 

0 

- 

0 

Al 

-.27132714 

KA 

,.02412880 

a2 

28175964 

a5 

0- 

a6 

0 

♦ 

B1 

.01462880 

Kb 

.01790107 

B2 

-1.9572863 

®5 

0 

B  6 

1.9260299 

Po 

P 

Cl 

.00874177 

KC 

.01933340 

C2 

.00896556 

c5 

0 

E 

P 

Al’ 

.  46473072 

.04286193 

A2' 

-.00620389 

V 

0 

m 

r 

Bl' 

-.02539885 

KfijRl2+  H2 

.03179911 

b2  ' 

.04309623 

b5’ 

0 

mm 

Al 

-.21275154 

ka 

.35205361 

a2 

. 46209967 

a5 

0 

a6 

0 

B1 

.01162773 

kb 

.26118594 

b2 

-3.2100413 

B5 

0 

b6 

3.1796150 

r0 

Cl 

.00685435 

Kc 

.28208436 

c2 

.01470395 

c5 

0 

- 

p 

Al’ 

.36440206  - 

%\K2  +  I12 

.62538141 

a2* 

-.01017468 

V 

0 

- 

r 

Bl’ 

-.01991605 

KbVB12  +  Ii2 

.46396553 

B2’ 

.07067985 

b5' 

0 

- 

0 

Al 

.3534235 

KA 

'.07815380 

a2 

-25.21345 

a5. 

0 

a6 

24.93682 

* 

B1  . 

-.01931556 

kb 

.05798158 

b2 

175.1489 

b5 

.6199628 

B6 

-175.1797 

C1C 

P 

Cl 

-.01138685 

KC  ' 

.06262090 

C2 

-.8022885 

c5 

.8679479 

- 

P 

Al’ 

-.60536104 

KaV^+Ii2 

; 13883429 

a2’ 

.55517272 

a5’ 

0 

- 

r 

Bl* 

.03308464 

Kb^  +  Ii2 

.1029990 

b2* 

-3.8565875 

b5* 

3.7886547 

-- 

0 

Al 

.07219731 

ka 

.1935925 

a2 

-16.45365 

a5 

0 

a6 

16.22009 

* 

B1 

-.00394581 

kb 

.1436248 

B2 

114.2976 

b5 

.4085555 

B6 

-114.1513 

Cnc 

P 

Cl 

-.00232607  ' 

KC 

,1551168 

c2 

-.5235526 

c5 

.3719777 

EE 

P 

Al' 

-.12366306 

KMtei2  +  Ii2 

.34390240 

a2’ 

,  . 36229131 

a5* 

0 

m 

r 

Bl’ 

.00675858 

Kb^.IiS 

.25513879 

B2* 

-2.5167086 

b5’ 

2.4967235 

M 

0 

Al 

.00235150 

ka 

.00311940 

a2 

.02525049 

a5 

0 

A6 

-.02886004 

+ 

B1 

-.00012851 

kb 

.00231425 

b2 

-.1754060 

b5 

0 

b6 

.1748510 

°*c 

P 

Cl 

-.00007576 

KC 

.00249943 

c2 

.00080347 

c5  . 

0 

- 

P 

Al’ 

-.00402776 

ka  K2  +  Ii2 

.00554138 

a2’ 

-.00055599 

V 

0 

- 

r 

Bl' 

.00022013 

Ksh2*1!2 

.00411110 

B2* 

.00386225 

b5’ 

0 

- 
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TABLE  II 

AMPLITUDE  COEFFICIENTS  OF  FREE  AND  TORCED  MOTIONS  OF  EXPERIMENTAL  SVEPT-WING  AIRPLANE  -  Concluded 


(b)  V  »  200  miles  per  hour. 


Initial  conditions 
and  forcing 
functions 

Variables 

Rolling- sub  s idence 
mode 

Oscillatory  mode 

Spiral  mode 

0 

Al 

0.01197011 

0.0l6l0600 

A  2 

0.4807800 

a5 

0 

A6 

0 

♦ 

®1 

-.00027^74 

.02347936 

b2 

-18.21875 

®5 

0 

b6 

18.21429 

0o 

P 

Cl 

-.0003213 

*c 

.02426556 

C2 

.0083529 

c5 

0 

- 

p 

Al’ 

-.02769022 

KA\jRi2  +  U2 

.03693084 

a2' 

-.00135277 

V 

0 

- 

*r 

Bl* 

.00049674 

KbK2+ii2 

.05285335 

b2’ 

.05126683 

b5' 

0 

-- 

■ 

0 

A1 

-.10082626 

Ka 

.13447276 

a2 

-.00492530 

a5 

0 

A6 

0 

* 

®1 

.00180874 

kb 

.19245072 

B2 

. 18665748 

b5 

0 

b6 

0 

& 

m 

.00270588 

Kc 

.19889500 

C2 

-.00008558 

c5 

0 

- 

P 

.23323073 

ka^ri2  +ii2 

.30270574 

A2’ 

.00001386 

V 

0 

- 

_ r _ 

-.00418398 

Kb^Ri2  +  IX2 

.43321690 

b2- 

-.00052525 

*>' 

0 

— 

D 

-.2093808 

% 

.00718425 

a2 

.2112909 

a5 

0 

A6 

0 

1 

E 

.00375619 

Kb 

.01028182 

b2 

-8.006691 

®5 

0 

B6 

7.995698 

Po 

P 

fm 

.005.61934 

Kc 

.01062610 

C2 

.00367088 

c5 

0 

- 

1 

.48433861 

KaVH12  +  I12 

.01617213 

A2’ 

-.00059456 

V 

■ 

- 

IBB 

Bl' 

-.00868879 

.02314494 

®2* 

.02253050 

B5' 

JM 

- 

A1 

-.07670895 

ka 

.1518581 

a2 

. 1869734 

a5 

i 

A6 

0 

(■■I 

B1 

.00137619 

Kb 

.2173312 

B2 

-7.085196 

b5 

1 

B6 

7.097827 

r0 

Cl 

.00205853 

KC 

.2246084 

c2 

.00324842 

1  1 

m 

. 

p 

Al' 

.17744277 

.34184086 

a2' 

-.00052614 

a5* 

1 

9 

1 

Bl’ 

-.00318338 

.48922468 

B2* 

.01993745 

b5' 

I 

■ 

A1 

. 4547069 

ka 

.03147098 

a2 

-365.6037 

a5 

0 

A6 

365.1805 

Bl 

-.00815719 

kb 

.04503932 

b2 

13855.46 

®5 

4. 457143 

B6 

-13855.50 

C>c 

Cl 

-.01220331 

Kc 

.04654752 

C2 

-6.351295 

c5 

6.400000 

H 

p 

Al’ 

-1.0518225 

Ka'JbI2-*'1!2 

.07084298 

a2  ' 

1.0286136 

a5* 

O' 

m 

r 

Bl* 

.01886912 

WRi2+Ii2 

.10138614 

®2* 

-38.981861 

b5* 

38.911304 

mm 

0 

A1 

.03526760 

kA 

.1276446 

a2 

-102.7051 

a5 

0 

a6 

102.5640 

♦ 

Bl 

-.00063272 

kb 

.1826785 

b2 

3892.267 

b5 

1.257143 

»6 

-3892.093 

C«c 

P 

Cl 

-.00094653 

Kc 

.1887956 

C2 

-1.784201 

c5 

1.600000 

1 

P 

V 

-.08158060 

kaJri2  +  ii2 

.28733535 

A2* 

.28895732 

a5* 

0 

fl 

r 

_JiL_ 

.00146360 

KB/r12  + !l2 

.41121949 

b2’ 

-10.950758 

b5* 

10.974984 

0 

A1 

.00140829 

ka 

.00193011 

a2 

.05655717 

a5 

0 

A6 

-.05882356 

♦ 

B1 

-.00002526 

kb 

.00276228 

b2 

-2.143383 

®5 

0 

®6 

2.142857 

c*= 

P 

C1 

-.00003779 

KC 

.00285477 

C2 

.00098269 

c5 

0 

m 

P 

Al' 

-.00325766 

Ka^i2  +  ii2 

.00434479 

A2 ' 

-.00015915 

a5' 

0 

r 

Bl' 

.00005844 

kb^i2  +  1i2 

.00621801 

V 

.00603139 

b5’ 

r,  radians  per  second,  p,  radians  per  second  j8,  radians  7^,  radians  0,  radians 


8.—  Response  of  experimental 
ero  time  of  a  constant  yawin 


